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Abstract. Einstein's field equations of gravitation are known to admit 
closed timelike curve (CTC) solutions. Deutsch approached the prob- 
lem from the quantum information point of view and proposed a self- 
consistency condition. In this work, the Deutsch equation is formulated 
as an eigenvalue problem. The disappearance of entanglement between 
two qubits in an Einstein-Podolsky-Rosen (EPR) state near a CTC is 
demonstrated. The method is utilized to analyze the discontinuous evo- 
lution of two chronology respecting (CR) qubits near a CTC. 

1 Introduction 

Godel pointed out the possibility of closed timelike curve (CTC) solu- 
tions for the Einstein's field equations of gravitation pQ. Several models 
and calculations gave evidence for existence of CTCs |2|3|4|5|6|7|8|9|10j . 
Because CTCs allow a time traveling particle to go back in time and inter- 
act with its own past, one can end up with paradoxes. Deutsch analyzed 
the problem from the quantum information point of view and proposed a 
self-consistency condition involving the density matrices of the chronology 
respecting (CR) and CTC qubits undergoing a unitary interaction (TT] . 
Deutsch showed that the self-consistency equation implies a non-unitary 
and nonlinear evolution for the CR components. Ralph introduced a toy 
model to show that unitarity can be recovered |12j . 

Consequences of existence of CTCs for quantum computation have 
been examined. Brun argued that using CTCs, composite numbers can 
be factorized efficiently with the help of a classical computer [13]. Bacon 
demonstrated that nonlinearity of Deutsch-like evolution can be used to 
efficiently solve computational problems by reducing NP-hard problems 
to P [14] . It has also been claimed that an observer with access to CTCs 
can perfectly distinguish nonorthogonal quantum states [15]. However, 
there are both counter and supporting arguments for this claim |16|17j . 

An alternative formulation of the quantum mechanics near CTCs is 
via post-selection [18119] . Lloyd et al, developed a model, which is phys- 



ically inequivalent to Deutsch's, based on combining quantum teleporta- 
tion with post-selection. Unlike Deutsch's mixed state generating CTCs, 
post-selected CTCs always send pure states to pure states, and hence, 
they do not create entropy. 

Pati et al, showed that in contrast to ordinary quantum mechanics 
where any mixed state can be purified by bringing an ancillary system, 
the states of quantum systems traveling in CTCs cannot be purified |20j. 
In other words, the CTC system interacting with the CR system cannot 
be viewed as a part of a larger CTC system in a pure entangled state 
in an enlarged Hilbert space. A practical consequence of this result is 
that in general it is not possible to simulate CTC qubits satisfying the 
Deutsch equation in laboratory, because arbitrary CTC states cannot be 
constructed in a consistent way. 

In this study, after introducing an eigenvalue formulation of the Deutsch 
equation, the method is used to demonstrate that the entanglement be- 
tween two qubits in Einstein-Podolsky-Rosen (EPR) state disappears 
when one of the two qubits interacts with a CTC qubit. The same for- 
mulation is utilized to examine the discontinuous evolution of CR qubit 
states, first pointed out by DeJonghe et al. [21]. 

2 The Deutsch Equation 

Deutsch's model involves a unitary interaction U of a CR system with 
another system that traverses a CTC. The equation is a result of the kine- 
matic self consistency condition, which implies that the density matrix 
of the CTC system after it interacts with the CR system is the same as 
that of the CTC system before the interaction. For each state of the CR 
register described by the density matrix pcr, CTC register is postulated 
to find a fixed point such that 



where <8> is tensor product and U* is the hermitian conjugate of U. Here, 
it is assumed that CR and CTC system are initially in a separable state. 
Deutsch showed that there is always a density matrix pctc which satisfies 
the above self-consistency condition by using the compactness of the space 
of density operators. In fact, there are in general many solutions, and 
one might need additional assumptions to decide which pctc to choose 
|11|14|21 22j. One of the immediate consequences of the equation is that 
the universe (i.e., the space of CR registers) may evolve from a pure state 




(1) 



to a mixed state. Nonlinear and noncontractive evolutions are among the 
other unusual effects |11|23|14] . 

Equation ([!]) implies that Pctc-, an d hence, its von Neumann entropy 

S (pctc) = ~Tr (pctc hi Pctc) (2) 

remains the same after CR and CTC systems interact. Therefore, the 
entanglement of pctc with the rest of the the enlarged CTC Hilbert 
space would remain the same for any local interaction if it were possible to 
distill pcTC, which runs counter to observation of Pati et al. [20J. However, 
for arbitrary local operations, the entropy does not remain constant but 
instead decreases. That is why, in spite of its simplicity, the demonstration 
of the Deutsch equation in laboratory is not a trivial task. One needs to 
find a two-particle scattering event, where the density matrix of one of 
the particles remains the same after the scattering. 



3 Eigenvalue Formulation 

Let CR and CTC registers be single qubits. For a pure state pen = 
without loss of generality pen can be assumed to be given by |0)(0| since 
U = U (V <£) I), where = V\0), is also unitary. For p C R = |0)(0| and 

\P21 P22 J 

inserting matrix elements of pen and Pctc, an d comparing the entries of 
its left and right hand sides, Eq. (JTJ) can be written as Mv = v where M 
is given in terms of the matrix elements it™ of unitary operator U, and 
their complex conjugates u*j, as 

/ UuU\ x + UziU% x UxiU*2 + ""31^32 u 12«ll + «32«31 n 12-"l2 + u 32«32 \ 
= uuu* 2l + u 31 u* 41 uuu* 22 + n 3 in| 2 ui 2 u* 2l + u 12 u\ 2 + u 32 u% 2 

^21^*1 + ^41""3l U\\U\ 2 + 1*41^2 U 2 2'«n + ^42^31 ^22""l2 + n 42«32 

\U 2 lU 21 + U4lU* 41 u 2 \u* 22 + u Al u\ 2 u 22 u* 21 + U42U41 u 22 u* 22 + u i2 u* 42 J 

and v = (pn, p\2, P21, P22) T ■ Here, v is normalized so as to satisfy p\\ + 
P22 = 1- Therefore, the Deutsch equation can be interpreted as a diago- 
nalization problem Mv = Xv, where one looks for A = 1 eigenvalues and 
vectors v whose entries define a valid density matrix. 

Although Deutsch gave a general proof for existence of a solution, it 
is worth rederiving the same result using the current formalism due to its 



simplicity. Using the unitarity of U, it is seen that the matrix M is of the 
form 



M = 



fa b b* c \ 

d e r 9 
d* f e* g* 

^1-a-b-b* 1-cj 



(5) 



Adding the last row of the matrix M — XI (I being the identity matrix) to 
its first row, one can immediately observe that 1 — A is always a factor of 
the characteristic equation det (M — XI). Therefore, there is at least one 
solution. If A = 1 is two-fold degenerate, any convex linear combination, 
av± + (1 — a) V2 with a € [0, 1], of the corresponding eigenvectors v\ and 
V2 is also a solution. As shall be seen below, there can be cases where 
i>2 does not correspond a valid density matrix, while its convex linear 
combinations with v\ are proper solutions. 

A practical way to construct M is to write it as M = A\\ ® A\ x + 
^21 ® ^21 where An and A21 are 2x2 matrices in 



and A*j denotes complex conjugation of Ay. 

Generalization to n CR qubits near a CTC qubit problem is straight- 
forward. In this case, 



where Aij are 2x2 matrices making up the 2 n+ x 2 n+ unitary matrix U 
as in Eq. ([6]). Even though the eigenvalue formalism has been presented 
for pure a pen, it can easily be generalized to mixed states. In this case 
M matrix is slightly more complicated, and it contains all four entries of 
PCR- 

4 Disappearance of Entanglement Near a CTC 

An interesting problem involving entanglement near a CTC is the be- 
havior of an EPR pair. In order to illustrate Deutsch's model, Bennett 
et al, considered half of a maximally entangled state, i.e., an EPR pair 
X (| 00) + |11}), and put it into a CTC [16J. Using both the single and 
multiple universe pictures, the authors showed that the joint state at 





M = ^ / A il ®A* : 



(7) 



i=i 



any time after the interaction between CR and CTC systems is a prod- 
uct state. In the eigenvalue formulation introduced in Sec. IV, the to- 
tal unitary interaction U involving three qubits (i.e., two CR qubits 
making up the EPR pair and one CTC qubit) can be written as U = 
(I 2 <g> E) (CNOT <g> I 2 ) (H <g> I 4 ) where H and CNOT denote the Hadamard 
and controlled-NOT operations, respectively. The EPR state is obtained 
by acting the operator {CNOT <g> I2) (H (g> I4) on the initial CR two qubit 
product state 1 00) , where 1% and I4 denote 2x2 and 4x4 unit matrices, 
respectively. Here, 



E 



/iooo\ 

00 10 
10 

yooo ly 



(8) 



denotes the exchange operation, where a qubit emerges from the CTC 
and half of the EPR pair is put into the CTC. 

Diagonalization of the corresponding M matrix gives the eigenvalues 
Ai = 1, and A2 = A3 = A4 = 0. Hence, the solution is unique, and the 
reduced density matrices for CR and CTC registers are I^/A and ^2/2, 
respectively. Therefore, two CR qubits are not entangled anymore. Since 
the reduced density matrix of the CR qubits is found to be mixed after 
the interaction, this is an expected result. It is known that according to 
the monogamy property of entanglement, if two qubits are maximally 
entangled, they cannot be entangled with a third qubit |26j . 



5 Discontinuous Evolutions Near CTCs 



DeJonghe et al, demonstrated that Deutsch's equation can lead to dis- 
continuities in the evolution of the CR systems [21] • The authors consider 
two CR qubits interacting with a CTC qubit via the unitary evolution 

U = |000}(100| + |100) <000| + |010}(011| + |011)(010| , 
|101)(110| + |H0)(101| + |001)(001| + I 111) <111 1 . ^ ' 

It is also assumed that the density matrix of CR qubits before the in- 
teraction is of the form pen = pi <8> P2- For three distinct initial states 
(Pcri Pcri Pgr) wm ch are infinitesimally close to each other, it has been 
shown that there is no choice of Pctc which is continuous in the vicinity 
of Pcr- Since the CR and CTC density matrices are obtained from the 
same matrix by partial trace operations, the discontinuity in the CTC 



state implies the behavior of the CR state. In the current formalism, 
starting again from CR state 1 00) , the three pure states Pcr-> Pcr-> Pcr 
correspond to the following M matrices 



M = 



/1-eOO e \ 
Oe 
e 
\ e 001-e/ 



(10) 



M 1 



/iooo\ 

00 00 
00 00 
\00 IJ 



(11) 



M c = 



(I 

v^O 7 

\0 







e 



1-eJ 



(12) 



where e G 
l-2e,\£ = -e,Xi 
for Pcr, and Af 



• A 3 



l,Af 



A 4 



[0, 1]. Eigenvalues for three cases are given by A^ 
'• 1 — e for initial state Pq R , Af = Af = 1, 

1,A£ = 1 - e,A 3 c = \<i = v^T^) for p c CR . 
Any eigenvector with eigenvalue 1 is a solution of the Deutsch equation. 
Degeneracy of Af and Af indicate that for initial state pQ R there are 
infinitely many solutions which can be obtained by taking convex linear 
combinations of the two degenerate states. CTC state solutions, all of 
which are independent of e, for the three cases are given by 



PCTC 




(13) 



PCTC 



/3 
1-/3 



for /3e [0,1], 



(14) 



c 

PCTC 




(15) 



Density matrices Pqj<q and p^tc can De obtained from Pq TC by choosing 
j3 = 1/2 and (3 = 1, respectively. As e — > 0, all three initial states approach 



each other, but Pq TC an< ^ Pctc being independent of e, are always differ- 
ent. Solving the Deutsch equation by a different method, DeJonghe et al. 
concluded that there is a discontinuity near p^ TC for the unitary interac- 
tion given above. Therefore, for finite e values, the eigenvalue formulation 
reproduces their result. However, as e — > 0, eigenvalues A^Af? — > 1, and 
hence the corresponding density matrices satisfy the Deutsch equation. 
At this limit, all three initial states lead to the same solution set, given 
by Eq. (fH|) . for pcrc- More explicitly, approximate solutions obtained 
by taking linear combinations of eigenvectors with eigenvalues Ai and A2 
for the initial states Pq R and p^ R can be written as 

P$tc = ( 1/2 ~ a 1/2 ° + a ) , for a £ [-1/2, 1/2], (16) 

and 

^TC=( 1 o 7 °),for 7 G[0,l]. (17) 

For both Pq R and Pc R , A 2 solution is an e independent vector V2 = 
(-1,0,0,1)' which gives a traceless 2x2 matrix. Even though this solu- 
tion cannot be associated with a density matrix, its linear combinations 
with v\ (eigenvector with eigenvalue Ai = 1), lead to proper vectors. So- 
lution set for all three cases at e — >■ is the same. For Pq R , the expression 
in Eq. (TH|) is a solution for all f3 € [0,1]. However, in case of pQ R and 
Pc R , the only e-independent and exact solutions are a = and 7 = 0, 
respectively. 

6 Conclusion 

The consistency condition proposed by Deutsch to avoid paradoxes near 
a CTC has been transformed to an eigenvalue equation. The proposed 
approach is a systematic method that solves the Deutsch equation. Two 
problems, namely an EPR pair near a CTC and the discontinuous evo- 
lution of the CR-CTC system, have been reexamined. Even though the 
eigenvalue formulation has been demonstrated for initially pure CR states, 
it can easily be generalized to the mixed case. 
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